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Our study of the energy distribution has shown that the strength of the first order transition in the four- 
dimensional compact U(l) lattice gauge theory decreases when the coupling A of the monopole term increases. 
The disappearance of the energy gap for sufficiently large values of A indicates that the transition ultimately 
becomes of second order. In our present investigation, based on a finite-size analysis, we show that already at 
A = 0.9 the critical exponent is characteristic of a second-order transition. Interestingly, this exponent turns out 
to be definitely different from that of the Gaussian case. We observe that the monopole density becomes constant 
in the second order region. In addition we find the rather surprising result that the phase transition persists up 
to very large values of A, where the transition moves to (large) negative /3. 



1. INTRODUCTION 

We investigate the compact U(l) lattice gauge 
theory in four dimensions with the Wilson action 
supplemented by a monopole term |^ 

5 = ^ (l-cose^,,,) + A^|Mp,,| (1.1) 

where Mp,^ = epativiQfj.j^.x+a -&iiu,x) /^t^ and the 
physical flux Q^u.x G [^""i is related to the 
plaquette angle 'dj^y.x £ (-47r,47r) by Q^iy^x = 
0Aii/,a; + 27rn^^ 3; g|. We use periodic boundary 
conditions. 

Our studies of the energy distribution |^,|| have 
shown that the strength of the first order transi- 
tion decreases with A so that ultimately the tran- 
sition becomes of second order. This observation 
has enabled us to develop an algorithm by which 
simulations have become possible also on big lat- 
tices 0. In this work we study the properties of 
the second order phase transition in detail. Be- 
cause the monopoles play a fundamental role in 
the dynamics of the phase transition, the addition 
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of a monopole term to the action appears partic- 
ularly well suited for exploring critical properties 
of the system appropriate for a continuum limit. 

2. PHASE TRANSITION LINE 

In Ref. iQl Barber and Shrock observed that 
there is a shift of the transition point if A is var- 
ied. The effects of a complete suppression of 
monopoles have also been studied H,^- In our 
previous investigations we have determined the 
location of the phase transition over a wide range 
of values in (/3, A) space. Now we wish to clarify 
what happens at very large A. For this purpose we 
determine the critical points at values of A sub- 
stantially larger than previously considered and 
allow (5 to take negative values. 

In order to keep the computational cost for this 
study within bearable limits we have used our 
topological characterization of the phases 1^,0. 
It is based on the fact that there is an infinite 
network of monopole current lines in the confin- 
ing phase and no such network in the Coulomb 
phase. On finite lattices "infinite" is to be de- 
fined in accordance with the boundary conditions 
1^ . For periodic boundary conditions "infinite" is 
equivalent to "topologically nontrivial in all direc- 
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tions" . Since the analysis of a single configuration 
(or of a few configurations, because of finite size 
effects, when close to the critical point) is already 
sufficient to identify the phase, this characteriza- 
tion permits to find the transition region rather 
quickly. 

It is to be remembered that on finite lattices 
different order parameters lead to slightly differ- 
ent critical /J. On an 8^ lattice the maximum of 
the specific heat and our topological order param- 
eter give values 1.0075(1) and 1.0074(2) for A = 0, 
and 0.3870(5) and 0.372(3) for A = 0.9, respec- 
tively. To determine the location of the maximum 
of the specific heat for larger A in an efficient way 
we first determine the critical (3 from the topolog- 
ical order parameter and then find the maximum 
of the specific heat in an easy second step. 

In Figure 1 we show our results for the location 
of the phase transition (3c , defined by the maxi- 
mum of the specific heat, for values of A ranging 
up to 1.3. It can be seen that the phase transition 
line continues to negative 13. 

Using the topological order parameter we could 
follow the line of phase transitions up to still much 
larger A: from A = 1.4 where I3c = —0.52(2) to 
A = 10 where f3c is approximately —1000. It is 
to be emphasized that the characteristic topolog- 
ical properties of the phases have been found to 
be fully present throughout the range of the in- 
vestigation. Thus we have the remarkable result 
that both phases are still present all the way up 
to very large A. 

We observe that finite size effects increase with 
A. This is indicated by the fact that the transition 
region becomes broader. The width of the peak 
of the specific heat increases by roughly a factor 
of 4 from A = to A = 0.9 and from A = 0.9 to 
A = 1.3. Similarly, from the topological order pa- 
rameter, we see an increase by a factor of approxi- 
mately 6 from A = to A = 0.9. A further indica- 
tion of larger finite size effects is that the height 
of the specific-heat maximum decreases with A. 
The decrease is roughly by a factor of 15 from 
A = to A = 0.9 and from A = 0.9 to A = 1.3. 
Thus one also sees that for very large A a precise 
determination of the location of the maximum of 
the specific heat becomes cumbersome. 



3. MONOPOLE NUMBER DENSITY 

At smaller A we have observed previously 
that the monopole number density in the 
Coulomb phase is roughly constant, while in the 
confining phase it decreases rapidly towards this 
constant with increasing A. Figure 2 gives the 
monopole number density along the transition 
line. It can be seen that starting at A = 0.9 the 
density becomes constant within the errors of the 
simulation. 

From the behavior exhibited in Figure 2 we 
conclude that the monopole number density be- 
comes constant, in the second order region. In- 
deed (as pointed out in Sect. 4) the critical be- 
havior characteristic of a second order transition 
properly occurs for A = 0.9 (while for A — 0.8 
there are still deviations from scaling). 



4. CRITICAL BEHAVIOR 

In order corroborate our observation that for 
large A the phase transition becomes of second or- 
der we have investigated the finite-size scaling be- 
havior of the maximum of the specific heat Cmax- 
It is expected to be 

Cmax L'' (4.2) 

if the phase transition is of first order and 

Cmax ~ (4.3) 

if it is of second order, where a is the critical 
exponent of the specific heat and v the critical 
exponent of the correlation length. 

In Figure 3 we present the simulation results for 
Cmax which we have obtained on lattices with L 
= 6, 8, 10, 12 for A = 0.9 at the corresponding 
values of [3c- The fit to these data gives 

- = 0.485(35) (4.4) 

Clearly this is quite far from 4 and thus the tran- 
sition not of first order. 
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By the hyperscaling relation a = 2 — dv we find 
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V = 0.446(5) 



(4.5) 



It is interesting to note that ( [4.5| ) is clearly dif- 
ferent from the value i of the Gaussian case. 



The value in Eq. (4.5) might be compared with 
V = 0.36(1) obtained recently in a sphere-like 
geometry. Further it should be remembered that 
values between 0.33 and 0.50 have been found, in 
different contexts, in works done in the eighties 

The critical /3 is expected to behave as 



(ic{L) = Pc{oo) + aL- 



(4.6) 



From this relation, using the value in Eq. (4.5) 
and our data for Pc{L) at A = 0.9, we get /3c(oo) 
= 0.4059(5) and a = -1.99(6). 

We have also performed simulations for A = 0.8 
on lattices with L = 6, 8, 10. It turns out that 
Cmax does not yet scale in this case. 

As a next step simulations at some value of A 
larger than 0.9 appear imp orta nt in order to see 
whether the value of Eq. (15) is universal. In 
any case our results suggest that a "line of fixed 
points along which a scale-invariant continuum 
limit will result" [|l2| occurs with the action (1.1) 
for sufficiently large values of A. 
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Figure 1. Location of phase transition Pc as func- 
tion of A for 8* (circles) and 16^ (crosses) lattices. 
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Figure 2. Monopole density as function of A. 
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Figure 3. (7max versus L for A = 0.9 at 13c- 



